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^ ! Abstract 
<N . 

The genuine quantum gravity effects can already be around us. It is likely that 
J> ■ the observed large-angular-scale anisotropies in the microwave background 

radiation are induced by cosmological perturbations of quantum-mechanical 
O . origin. Such perturbations are placed in squeezed vacuum quantum states 

and, hence, are characterized by large variances of their amplitude. The 
\ statistical properties of the anisotropies should reflect the underlying statis- 

^ ■ tics of the squeezed vacuum quantum states. In this paper, the theoretical 

O . variances for the temperature angular correlation function are described in 

I . detail. It is shown that they are indeed large and must be present in the ob- 

servational data, if the anisotropies are truly caused by the perturbations of 
quantum-mechanical origin. Unfortunately, these large theoretical statistical 
uncertainties will make the extraction of cosmological information from the 
measured anisotropies a much more difficult problem than we wanted it to 
be. This contribution to the Proceedings is largely based on references [42,8]. 
The Appendix contains an analysis of the "standard" inflationary formula for 
density perturbations. 
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I. INTRODUCTION 



In the context of cosmology, the quantum gravitational physics is usually understood as 
the early Universe physics at the Planck scale. It is assumed that the gravitational field 
is fully quantized and the Universe itself is, in a sense, quantized. (For a comprehensive 
recent review of the structural issues in quantum gravity see [1]). However, there exists also 
another meaning in which we can speak of quantum gravity effects. To clarify the difference, 
let us consider an analogy from condensed matter physics. Imagine a crystal and various 
quantized excitations in it: phonons, rotons, excitons, etc.. The creation and annihilation 
operators that the condensed matter theorists write do not create and annihilate the crystal, 
they create and annihilate excitations in the crystal. The excitations should not necessarily 
be linear, one can take account of the higher-order corrections too. The theory of the 
crystal excitations is a fully legitimate quantum theory which does not attempt to quantize 
"everything in sight". Similarly, in our study below, we do not write the creation and 
annihilation operators that create and destroy universes. Our operators create and destroy 
perturbations in the Universe. Nevertheless, the effects that we are studying are genuine 
quantum gravity effects in the sense that they inherently contain all the three fundamental 
constants. The gravitational constant and velocity of light enter because we deal with a 
gravitational field (its energy-momentum tensor), the Planck constant enters because we 
normalize the vacuum energy of the field to have "a half of the quantum" in each mode. All 
three fundamental constants combine in the Planck length Ipi, and this quantity naturally 
appears as the coefficient in the most of our formulas. The sending of the Ipi to zero would 
eliminate the entire expression. 

Now, why at all do we think that the anisotropies in the microwave background radiation 
may have something to do with gravitational quantum physics? Why, in the first place, is 
there such a considerable interest to the measured large-angular-scale anisotropies [2]? After 
all, we have always knew that the observed part of the Universe is homogeneous and isotropic 
only approximately, on average, and we are certainly aware of large deviations on smaller 
scales. 

The point is that the photons of the microwave background have been traveling to us 
for almost all duration of that characteristic time which we call the age of the Universe 
and which is determined by the value of the measurable quantity — the Hubble parame- 
ter. The anisotropies on largest angular scales are specific in that they are prodiiced by 
cosmological perturbations on largest spatial scales. In our context, the largest spatial scale 
means definitely larger than any directly studied distance, and comparable or larger than 
the characteristic cosmological distance — the Hubble radius — associated with the Hubble 
parameter. In terms of the crystal analogy, and assuming that the size of the crystal varies 
at some time scale, we are dealing with excitations whose wavelengths are comparable and 
longer than the light travel spatial scale associated with the time variation. 

The next question is whether we will be attempting to explain the origin and nature of 
such long-wavelength cosmological perturbations or we will be happy to simply accept their 
existence. There is no logical contradiction in the second position. One can study the per- 
turbations with whatever accessible observational accuracy one has, and then extrapolate 
back in time their evolution according to classical dynamical laws. Typically, we will end up 
with a very anisotropic and inhomogeneous universe at some very early times. As an expla- 
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nation, we will have nothing more to say except that this was the Universe that was given 
to us for our future life and study. A more appealing position, at least aesthetically, is to try 
and find out a universal and reliable mechanism which could be capable of generating the 
required perurbations in the originally homogeneous and isotropic Universe. The problem 
is not to increase the wavelengths of the perturbations up to the size of the present Hub- 
ble radius (in the expanding Universe, the wavelength is always growing with time) but to 
generate them. In principle, one can imagine that such a mechanism may have operated in 
the relatively recent Universe and may have not be related to quantum gravity. In practice, 
taking into account everything what we know about cosmology and the necessity to produce 
perturbations with extremely long wavelengths, it is difficult to suggest such a mechanism, 
especially if we do not want to make many additional hypotheses. It appears that we can 
only rely on the quantum processes in the very early Universe. The parametric interaction 
of quantized cosmological perturbations with strong variable gravitational field of the very 
early Universe provides us with such a possibility. The perturbations are generated as a 
result of amplification of their zero-point quantum oscillations. Returning to the laboratory 
physics analogy, one can recall that a cavity filled with a dielectric medium and initially free 
of electromagnetic radiation will eventually contain the radiation if the parameters of the 
dielectric medium vary properly in time. The quantum-mechanical (parametric) mechanism 
of generating cosmological perturbations relies only on the validity of the general relativity 
and the basic principles of quantum field theory. The observational consequences of this 
phenomenon we will study below. 

The line of reasoning in this study can be summarized as follows. 

We see the anisotropics in the microwave background radiation at the largest angular 
scales [2]. Observers convincingly argue that this is a genuine cosmological effect. 

If the large- angular-scale anisotropy in the microwave background is really produced by 
cosmological perturbations (density perturbations, rotational perturbations, gravitational 
waves), then their today's wavelengths are of the order and longer than today's Hubble 
radius In- Strictly speaking, all wavelengths give contributions to the anisotropy at every 
given angular scale. But if the spectrum of the perturbations is not excessively "red" or 
"blue" , the dominant contribution is provided by wavelengths indicated above. For instance, 
the major contribution to the quadrupole anisotropy is provided by wavelengths somewhat 
longer than Ih- 

In the expanding Universe, the wavelengths of perturbations increase in proportion to 
the cosmological scale factor. The wavelengths that are longer than some length scale 
today have always been longer than that scale in the past. Moreover, the wavelengths of the 
perturbations of our interest are much longer than the Hubble radius defined at the previous 
times, when one goes back in time up to the era of primordial nucleosynthesis — the earliest 
era of which we have observational data. It is hard to imagine (although it does not seem 
to be logically impossible) that cosmological perturbations of our interest, with such long 
wavelengths, could have been generated by local physical processes during the interval of 
time between the era of primordial nucleosynthesis and now. We are bound to conclude that 
these perturbations were generated in the very early Universe, before the era of primordial 
nucleosynthesis. There is still 80 orders of magnitude, in terms of energy density, to go from 
the era of primordial nucleosynthesis to the Planck era; a lot of things could have happened 
in between. 



3 



The law of evolution of the very early Universe is not known, but it is likely that it 
could have been significantly different from the law of expansion of the radiation-dominated 
Universe. If so, some amount of cosmological perturbations must have been generated 
quantum-mechanically, as a result of parametric interaction of the quantized perturbations 
with strong variable gravitational field of the very early Universe. Gravitational waves have 
been generated inevitably, while density and rotational perturbations — if we were lucky; 
see [3] and references therein. (If the cosmological scale factor has always been the one 
of the radiation-dominated Universe, we must stop here, because the parametric coupling 
vanishes in this case, and cosmological perturbations cannot be amplified classically and 
cannot be generated quantum-mechanically.) The amount and spectrum of the generated 
perturbations depend on the law of evolution of the very early Universe (the strength and 
variability of the gravitational pump field), and this is how we can learn about what was 
going on there. In particular, the law of evolution of the very early Universe could have 
been of infiationary type. 

If the cosmological perturbations were generated quantum-mechanically, they should be 
placed in the squeezed vacuum quantum states [4] (for an introduction to squeezed states see, 
for example, Ref. [5,6] and the pioneering works quoted there). Squeezing of cosmological 
perturbations might have degraded by now at short wavelengths but should survive at long 
wavelengths, especially in the case of gravitational waves. 

The squeezed vacuum quantum states can only be squeezed in the variances of phase 
which unavoidably means the increased variances in amplitude. The statistical properties of 
the squeezed vacuum quantum states are significantly different from the statistical properties 
of the "most classical" quantum states — coherent states. This is well illustrated by the fact 
that the variance of the number of quanta in a strongly squeezed vacuum quantum state 
is much larger than the variance of the number of quanta in the coherent state with the 
same mean number of quanta (N), (N) ^ 1. For a squeezed vacuum state the variance is 
(N^) - {Nf = 2{N){{N) + 1) > (iV), while for a coherent state it is (N^) - {Nf = (N). In 
cosmology, the mean number (A^) is a characteristic of the expected mean square amplitude 
of cosmological perturbations, while the variance (A^^) — (iV)^ is a characteristic of theoretical 
uncertainties in the amplitude. These two characteristics are independent properties of a 
quantum state or a stochastic process. Theoretical models may agree on (N) and disagree 
on variance or agree on variance and disagree on (A^). 

The statistical properties of squeezed cosmological perturbations will inevitably be re- 
fiected in statistical properties of the microwave background anisotropics caused by them. 
Squeezing is a phase-sensitive phenomenon, and to fully extract its properties the quantum 
optics experimenters use the phase-sensitive detecting techniques based on a local oscillator. 
In cosmology, we are very far from being able to build a local oscillator, except of maybe, 
in the distant future, for short gravitational waves. Besides, in our study of the microwave 
background anisotropics, we are interested in so long-wavelength perturbations that it would 
take billions of years to wait for seeing the time dependent oscillations of variances in the 
quadrature components of the perturbation field. On the other hand, the amount of cosmo- 
logical squeezing is enormously greater than what is achieved in quantum optics laboratory 
experiments. In cosmology, we can only rely on the phase-insensitive, direct detection. One 
can expect that the underlying large variances of the amplitude of cosmological perturba- 
tions should result in large statistical deviations from the mean values for the microwave 
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anisotropies. 

A detailed study and proof of this statement is the purpose of this work. 

At this point it is necessarry to comment on the possible numerical values of {N) (and, 
hence, the amplitude) for cosmological perturbations of different nature which can be gener- 
ated quantum-mechanically in one and the same cosmological model. A consistent quantum 
theory provides us, of course, with both, the mean value of and its variance. According 
to the calculations of Ref . [3] , the contribution of quantum-mechanically generated gravita- 
tional waves to the large-angular-scale anisotropy is somewhat greater (even in the limit of 
the de Sitter expansion) than the contribution of quantum-mechanically generated density 
perturbations. It is argued in [7,8] that the "standard" inflationary formula for density per- 
turbations, which requires in this limit an arbitrarily large excess of density perturbations 
over gravitational waves, is based on errors. Wc additionally discuss this issue in some de- 
tail in the Appendix of this paper. However, the major emphasis of this paper is on the 
statistical properties of cosmological perturbations which are determined by their common 
origin — quantum mechanics and squeezing. These properies are related to the variance 
of A^ rather than to its mean value. Our discussion will be equally well applicable to the 
perturbations of any nature (density perturbations, rotational perturbations, gravitational 
waves) if they have the same origin. 

II. THE GENERAL EQUATIONS FOR QUANTIZED COSMOLOGICAL 

PERTURBATIONS 

Here we will briefly summarize some basic information about quantized cosmological 
perturbations (see [3,7] and references therein). The squeezed field operator derived in 
this Section is a basic mathematical construction for our further discussion of statistical 
properties. 

The metric of the homogeneous isotropic universe can be written in the form 



where is the spatial metric. For reasons of simplicity, we will be considering only spatially 
flat universes, that is jij = Sij. 

Following Lifshitz, it is convenient to write the perturbed metric in the form 



where hij are functions of r^-time and spatial coordinates. By writing the perturbed metric 
in this form we do not lose anything in the physical content of the problem, but we gain 
considerably in the mathematical tractability of the perturbed Einstein equations. The one 
who is interested in solving equations will certainly be interested in their simpler form. 
Those who prefer "gauge-invariant formalisms" are welcome to take the found solution and 
compute with its help whichever gauge-invariant quantity they like. These quantities, being 
gauge-invariant, have the same values in all gauges. 

The components h^j of the perturbed gravitational fleld can be classifled in terms of 
scalar, vector, and tensor eigenfunctions of the Laplace differential operator. The compo- 
nents of the perturbed energy-momentum tensor can also be classifled in the same manner. 



ds"^ — —a^{r]){dr]^ — ■fij dx^ dx^) 



(1) 



ds^ = -a'^{'q)[drf - {5ij + hij)dx'' dx^] 
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After that, the hnearized Einstein equations reduce to a set of ordinary differential equa- 
tions, separately for scalar (density perturbations), vector (rotational perturbations), and 
tensor (gravitational waves) parts. 

The number of independent unknown functions of time that can potentially be present 
(on grounds of the classification scheme) in the perturbed Einstein equations is always 
greater than the number of independent equations. It is 6 functions and 4 equations for 
density perturbations, 3 functions and 2 equations for rotational perturbations, and 2 func- 
tions and 1 equation for gravitational waves. In order to make the system of equations closed, 
it is necessary to say something about the perturbed components of the energy-momentum 
tensor or to specify from the very beginning the form of the energy-momentum tensor. The 
popular choices are perfect fluids and scalar fields. Even for gravitational waves, it is not a 
totally trivial question what their definition is (see, for example, Ref. [9]). However, after 
everything is being set, and as soon as the scale factor a{r]) (the background solution) is 
known, the general solution to the perturbed equations can be found. In practice, exact so- 
lutions are being found piecewise, at the intervals of evolution where the energy-momentum 
tensor has simple prescribed forms. 

We can now write the quantum-mechanical operator for the perturbations of the gravi- 
tational field hij in the following universal form: 



C 
a{r]) 



1 r°° ^ s 1 



c„(r7)e-- + c:(r;)e- 



(3) 



We will start the explanation of Eq. (3) from the polarization tensors Pij. Let us intro- 
duce, in addition to the unit wave- vector n/n, two more unit vectors /j, rrii, orthogonal to 
each other and to n: 



Hi 

n 



(sin 6 cos 0, sin 6 sin 0, cos 6) , /j 
±(cos cos 0, cos 6 sin 0, — sin 6) , 



(sin 0, — cos 0, 0) , 



(4) 



+ for ^ < f , - for e>l. 

The two independent polarization tensors, s — 1,2, for each class of perturbations, can 
be written as follows. For gravitational waves: 

1 2 

Pij{n) = - miinj) , Pii(n) = {liVLj + Ijiiii) . 
For rotational perturbations: 

1 1 2 1 

Pij{n) = -{huj + IjUi) , Pij{^) = -{rriinj + rrijUi) . 

Th Th 

For density perturbations: 



n) 




^J 1 



-V3^ 



The polarization tensors of each class satisfy the conditions PijP^^ — 2Sss', Pij{—''^) = Piji'^)- 
In practical handling of the density perturbations it proves convenient to use sometimes, in 
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1 

addition to the scalar polarization component Pij, the longitudinal-longitudinal component 

2 

(proportional to riiTij) instead of Pij. The explicit functional dependence of the polarization 
tensors is needed for the calculation of various angular correlation functions. 

The evolution of the creation and annihilation operators Cn(77), 0^(77), for each class of 
perturbations and for each polarization state, is defined by the Heisenberg equations of 
motion: 

—j^ = -i[cr,[ri),H] , = -t[ci{r]),H] . (5) 

The dynamical content of the problem is determined by the Hamiltonian H. Its form depends 
on the class of perturbations and additional assumptions about the energy-momentum tensor 
which we have to make, as was discussed above. 

Under the simplest assumptions about gravitational waves (waves interact only with the 
background gravitational field, there is no anisotropic material sources) the Hamiltonian for 
each polarization component takes on the form 

H = nclcn + nc^nC-n + 2(7(77) cJ^cLn + 2(7*(?7)c„c_„ (6) 

where the coupling function a{r]) is a{r)) = | ^■ 

For rotational perturbations, assuming that the primeval matter is capable of supporting 
torque oscillations, assuming that the oscillations are minimally coupled to gravity, and 
assuming that the torsional velocity of sound is equal to the velocity of light, the Hamiltonian 
for each polarization component reduces to exactly the same form (6) with the same coupling 
function a{f]). 

For density perturbations, we consider specifically a minimally coupled scalar field with 
arbitrary scalar field potential as a model for matter in the very early Universe, and perfect 
fluids at the later eras. The quantization is based on the scalar polarization component (the 
function of time responsible for another polarization state is not independent). There is 
only one independent sort of creation and annihilation operators in this case. The operators 

Cn(?7), c^{f]) are expressible in terms of the operators d^irj), dl^{ri) for which the Hamiltonian 
has again the same form (6) but with the coupling function a{7]) = f-^^^-, where 



For density perturbations, it is the operators dn{r]), dl^{rj) that participate in Eqs. (5), (6). 

Now, let us turn to the constant C in Eq. (3). Its value is determined by the normalization 
of the field of each class to the "half of the quantum in each mode" . Under the assumptions 
fisted above, one derives C = VIStt Ipi for gravitational waves, C = \/32ttIpi for rotational 
perturbations, and C = -\/247r Ipi for density perturbations, where Ipi is the Planck length, 

ipi = {Gh/^Y'\ 

The form of the Hamiltonian (6) dictates the form of the solution (Bogoliubov transfor- 
mation) to Eq. (5): 

Cn(^) = Un{r])Cn{0) + Vn{ri)c^_^{{)) 

cU^) =<(^)cL(0)+<(r7)c_n(0) (7) 
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where Cn(0), cJj(O) are the initial values of the operators taken long before the interac- 
tion with the pump field became important (a(ri)/n 0) and which define the vacuum 
state Cn(0)|0) = 0. The complex functions Un{r]), Vn{r]) obey coupled first-order differential 
equations following from Eq. (4) and satisfy the condition — \vn\'^ — 1 which guaran- 
tees that the commutator relationship [cn(0), cJji(O)] = S^{n — m) is satisfied at all times, 
[cn(77), cjjj(77)] = 5^{n — m). If one introduces the function fin{v) = Unifj) + v^{r]), one recov- 
ers from the equations for Vn{r)) the classical equations of motion. For gravitational 



waves: 



n — 



(8) 



For rotational perturbations: 



n 



a 
a 



(9) 



where Vt is the torsional velocity of sound which we assumed above to be c. For the scalar 
field density perturbations: 



n 



(10) 



If the pump field is such that the 7 function is independent of time, Eq. (10) reduces to 
exactly the same form as Eq. (8) for gravitational waves. 

In the Schrodinger picture, the initial vacuum quantum state |0n) |0-n) evolves into a 
two-mode squeezed vacuum quantum state. In our problem, each of the two-mode squeezed 
vacuum quantum states is a product of two identical one-mode squeezed vacuum quantum 
states which correspond to the decomposition of the real field hij over real spatial harmonics 
sinnx and cosnx. In the Heisenberg picture, the initial vacuum quantum state does not 
evolve in time and is the same now. 

By using Eq. (7) one can present the field (3) in the form 



1 ^ s 1 



.(0) 



where the functions hniv) hniv) = '^^b^niv) + '^niv)]- For gravitational waves and 

s s s s 

rotational perturbations, the functions hn are simply hn = f^n/o- where are solutions to 
Eqs. (8), (9) with appropriate initial conditions. For density perturbations, the functions 

s 

hn are derivable from solutions to Eq. (10) in accord with the relationship between c and d 

1 2 I'l" 2I" 

operators. Besides, for density perturbations, we should regard c^{0) = c^{0), Cj^(O) = c„(0) 

in Eq. (11). In all cases, for a given cosmological model, that is for a model in which the 

s 

scale factor 0(77) is known from the very early times and up to now, the functions hn can be 
found from the classical equations of motion with appropriate initial conditions. 

It follows from Eq. (11) that the mean quantum-mechanical value of the field h^j is 
zero at every spatial point and at every moment of time, (0|/iij|0) = 0. One can also 
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calculate variances of the field, that is the expectation values of its qiiadratic combinations. 
One useful quantity is hijh^K By manipulating with the product of two expressions (11), 
using the summation properties of the polarization tensors, and remembering that the only 
nonvanishing correlation function is 

(0|c„(0)cl,(0)|0) = 5,,,5=^(n-n') , 

one can derive the formula 

(0|/i,,(r?,x)/i'^(7?,x)|0) = — nY^hn{v)?dn . (12) 

Equation (12) shows that the variance is independent of the spatial point x but does depend 
on time. 

The expression under the integral in formulas such as Eq. (12) is usually called the power 
spectrum (in this case, it is the power spectrum of the quantity hijh^^): 

Pin) ^ ^,nj2\Uv)\' ■ (13) 

In cosmology, it is common to use the power spectrum defined in terms of the logarithmic 
frequency interval, that is the function 

Pz{n) = ^yj:\Uv)\' (14) 

{Z from Zeldovich). We are mostly interested in the power spectrum of cosmological per- 
turbations in the present Universe, at the matter-dominated stage. This spectrum is never 
smooth as a function of frequency (wave-number) n. Squeezing and associated standing 
wave pattern of the field make the spectrum an oscillating function of n for each moment 
of time. In their turn, the oscillations in the power spectrum will produce oscillations in 
the distribution of the higher-order multipoles of the angular correlation function for the 
temperature anisotropies. However, the spectrum is smooth for sufficiently long waves. At 
a given moment of time, it applies to all perturbations whose wavelengths are of the order 
and longer than the Hubble radius defined at that time. Moreover, the smooth part of the 
spectrum is power-law dependent on n if the scale factor a{r]) of the very early Universe (the 
pump field) was power-law dependent on 77-time. 

Let us assume that the scale factor at the initial stage of expansion was 

a{rj) = /o|r7|^+^ (15) 

where lo and j3 are constants. If the evolution is governed by a scalar field, the Einstein 
equations require the constant (3 to he (3 < —2. The value f3 = —2 corresponds to the 
de Sitter expansion. At later times, the scale factor changed to the laws of the radiation- 

s 

dominated and matter-dominated universes. From solutions for h{i]) traced up to the matter- 
dominated stage, one can find 
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iLivW - r2'''^^' and Pz{n) ^-^n'^f'^'^ 



It is convenient to introduce the characteristic amplitude h{n) of the metric perturbations 
defining this amplitude as the standard deviation (square root of variance) of the perturbed 

gravitational field per logarithmic frequency interval. In the long-wavelength limit under 
discussion, this quantity is universally expressed (both, for gravitational waves and density 
perturbations) by the formula [21]: 

Mn)~y^n^+' . (16) 

I'D 

Note that the functional form of h{n) is the same for gravitational waves and density per- 
turbations, the difference is in the numerical coefficient (omitted in this discussion) which 
is somewhat in favor of gravitational waves [3,7]. The numerical level of h{n) is mainly 
controlled by the constant lo- 

The spectra of other quantities can be found in the same manner. For instance, in case 
of density perturbations, one can derive the spectrum of perturbations in the matter density 
Sp/p. Since the relationship between Sp/p and the metric perturbations involves the factor 
(nr))'^ and, hence, involves two extra powers of n, ^(n) ~ n'^h{n), one finds 

I V f°° DP/ 



p p Jo n 



where P|(n) ~ and 



Sp, . Ipi ^+4 , s 

— (n) ~ —n^^ . (17) 
P h 

It follows from Eq. (16) that h{ri) is independent oi n ii P — —2. This independence 
corresponds to the original Zeldovich's definition of the "flat" spectrum: all waves enter 
the Hubble radius with the same amplitude. If the gravitational field perturbations h{n), 
regardless of their wavelength, have equal amplitudes upon entering the Hubble radius, the 
matter density perturbations ^{n) do also have equal amplitudes (the extra factor [nr])'^ is 
of the order of 1 at the time when a given wave n enters the Hubble radius) . For models of 
the very early Universe governed by a scalar field, the spectral index P + 2 in Eq. (16) can 
never be positive. 

Formula (16) and the associated formula (17) should be compared with the "standard" 
inflationary formula which requires that the amplitudes of density perturbations taken at 
the time of entering the Hubble radius should go to inflnity in the limit of the de Sitter 
inflation, /3 — > —2. It should also be noted that a "disgusting convention" (the term is 
borrowed from Ref. [10]) is often being used according to which one and the same Harrison- 
Zeldovich spectrum is described by the spectral index Ut = for gravitational waves and 
by the spectral index = 1 for density perturbations. Of course, there is no need in this 
convention. In both cases, the metric perturbations with the Harrison-Zeldovich spectrum 
are described by the same spectral index (zero), see Eq. (16). 

We will finish this section with a short discussion of coherent states. There is no natural 
mechanism for the generation of cosmological perturbations in coherent states, but if there 



10 



were one it would be reflected in many parts of the theory. The interaction part of the 
Hamiltonian (6) would be linear (not quadratic) in the creation and annihilation operators. 
The analogue of Eq. (7) would read 



where the complex function an{f]) is determined by the coupling function in the Hamiltonian. 
On the position — momentum diagram, the evolution (18) of the field operators corresponds 
to displacing the vacuum state without squeezing whereas the evolution (7) corresponds 
to squeezing the vacuum state without displacing. In terminology of mechanics, coherent 
states arc produced by a force acting on the oscillator whereas squeezed vacuum states 
are produced by a parametric influence. In coherent states, the mean value of the field is 
not zero. The correlation functions (at least, for some quantities) would also be different 
from the squeezed state case. In cosmological context, this would eventually be reflected 
in the differing statistical properties of perturbations and induced microwave background 
anisotropics. 

The calculations of the next Section are based in an essential way on the field operator 
(11) for the squeezed vacuum perturbations. 

III. QUANTUM-MECHANICAL EXPECTATION VALUES FOR THE 
MICROWAVE BACKGROUND ANISOTROPIES 

The microwave background anisotropics are a subject of intense study [22]. 

In absence of cosmological perturbations, the temperature of the microwave background 
radiation seen in all directions on the sky would be the same, T. Let us denote a direction 
on the sky by a unit vector e. The presence of cosmological perturbations makes the tem- 
perature seen in the direction e differing from T. The temperature perturbation produced 
by density perturbations or gravitational waves can be described by the formula [11]: 



where dhij/drj is taken along the integration path = e*w, f) = tJr — w, from the event of 
reception w — to the event of emission w — Wi — rjn — tje- The formula for rotational 
perturbations is more complicated than (19) [11], and we will leave rotational perturbations 
aside. 



c„(77)=e-^"^c„(0)+a„(?7) 



(18) 




(19) 



For quantized cosmological perturbations, the 
erator. Using Eq. (11) we can write this operator as 



(e) becomes a quantum-mechanical op- 



ST 




X 



S„(0)/J^)e^-- + cl(0)/l(«;)e-^' 



(20) 



where 
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V=Vr-w 



Having defined the observable ^(e) and knowing the quantum state |0) we can compute 
various quantum-mechanical expectation values. In the laboratory quantum mechanics, 
the verification of theoretical predictions expressed in terms of the expectation values would 
require experiments on many identical systems. An immediate generalization of this principle 
to cosmology would require speculations about outcomes of experiments performed in "many 
identical universes". Without having access to "many universes" we can only rely on the 
mean (expected) values of the observables and on the probability distribution functions as 
indicators of what is likely or not to be observed in our own single Universe. We will return 
to this point in Sec. IV. 

The expected value of the temperature perturbation to be observed in every fixed direc- 
tion on the sky is zero: 

(0|f^(e)|0)=0 . 

One particular measured temperature map is the result of the measurement performed 
over one particular realization of the random process describing cosmological perturbations 
of quantum-mechanical origin. For this realization, the temperature perturbations may, 
should, and in fact are, present. Many measurements will not help (except of reducing the 
instrumental noises) in the sense that they all should give identical results, because the 
timescalc of the perturbations under discussion is so enormously larger than an interval of 
time between the experiments. If the COBE's map is correct, we will have to live with this 
map practically forever. 

Let us now compute the expected angular correlation function for the temperature per- 
turbations seen in two given directions on the sky, ei and 62. This correlation function is 
defined as the mean value for the product of ^(ei) and ^"(^2): 

^(ei,e2) = (0|^(e0^(e2)|0) . (21) 
By manipulating with the product of two expressions (20) one can derive the formula 

1 1 rwi rwi rco 

K(ei,e2) ^ -C^T^ dw dw d'^ne'"'^^''"-^^''^ 
4 (27r)3 Jo Jo J-00 

X E (P^,(n)elei) [kA^)e^,ei) fMfliw) . (22) 



The next step is the formidable task of taking the integrals over angular variables in 3- 
dimensional wave-vector n space. However, it can be done (see Ref. [12] for gravitational 
waves, Ref. [13] for rotational perturbations, and Ref. [3] for density perturbations). The 
final expression reduces, without making any additional assumptions whatsoever, to the 
form 

00 

X(ei, 62) = K{S) = III E PiicosS) . (23) 
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Wc sec that the correlation function depends only on the angle S between the directions ei, 62 
not directions themselves. The coefficient /p; is taken from C^, other numerical coefficients 

s 

are included in Ki. The quantities Ki involve the integration of fn{w) over the parameter w 
and the remaining integration over the wave-numbers n. The numerical values of Ki depend 
on a chosen sort of cosmological perturbations and a chosen cosmological model; so far, the 
formula (23) is totally general. Pi{cos 5) are the Legendre polynomials. The lowest multipole 
^min follows automatically from the theory and it turns out to be, not surprisingly, /min = 
for density perturbations, /min = 2 for gravitational waves (and /min = 1 for rotational 
perturbations). For the separation angle 5 = 0, Eq. (23) reduces to the variance of ^(e), 
that is 



(0|f^(e)f^(e)|0) = X(0) 



''PI 



(24) 



Formula (23) gives the expected value of the observable ^ (ei) ^ (62) . If the experimenter 
measured this observable in "many universes" and averaged the measured numbers, he/she 
would get the result (23). Moreover, formula (23) says that if the experimenter made the 
measurements at any other pair of directions, but with the same separation angle 5, he/she 
would again get, after the averaging over "many universes", the result (23). Without having 
access to "many universes", we can ask what is the theoretical standard deviation of the 
quantity ^(61)^(62). (In practice, for deriving K{5) we need a kind of ergodic hypothesis 
allowing us to replace the averaging over "universes" by the averaging over pixels on a single 
map.) The variance V(ei, 62) of this quantity is, by definition. 



nei,e2) = (0|^(e0^(e2)^(e0^(e2)|0) - 



(0|f^(ei)^(e2)|0) 



(25) 



The standard deviation is the square root of this number. 

The calculation of ^(ei, 62) requires us to deal with the product of four expressions (20). 
However, the mean values of the products of four creation and annihilation operators are 
easy to handle. One can show that 1^(61,62) depends only on the separation angle 5 and 



V{e^,e2)=V{5) 



, 5T , , 6T , , , , 
I ^(ei) — (62)10) 



+ 



, , 5T , , 5T , , , , 
(0| — (e) — (e)|0) 



that is 

The standard deviation for the observable ^"(61)^(62) is 



(26) 



(27) 



(28) 



In a similar fashion one can derive the higher order correlation functions for two directions 
ei, 62 and the correlation functions for larger number of directions, but we will not need 
this information. 

In the limit 5 = Eqs. (25), (26) say that 
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2 

(29) 

The familiar factor 3 relating the fourth-order moment with the square of the second-order 
moment (given that the first-order moment is equal to zero) is the reflection of the underlying 
Gaussian nature of the squeezed vacuum wavefunctions associated with the Hamiltonian (6). 

By examining Eq. (28) one can conclude that for each separation angle 5 the standard 
deviation of the angular correlation function is very big. Even at those separation angles 
at which K{6) vanishes, the standard deviation is as big as the variance for ^(e) itself. 
However, the value of the standard deviation for a given variable is not very informative 
per se, as long as the probability density function for this variable is not known. If the 
probability density function (p.d.f.) were normal, we could say that the probability to find 
a result outside of la interval is 32%. Without knowing the p.d.f. we could resort to the 
Chebyshev inequality, but it would only tell us that this probability is less than 1. To 
get more information about possible deviation of the angular correlation function from its 
mean values we will consider in Sec. IV a classical random model which will reproduce the 
expectation vahics calculated above and will allow us to construct the p.d.f. for the variable 
^(ei) ^"(62)- On the other hand, the quantum-mechanical calculations of this Section will 
shed light on the classical model. As is known, "quantum mechanics helps us understand 
classical mechanics", see on this subject a paper of Zeldovich signed by the pseudonym 
Paradoksov [14]. 



(0| 



'6T 



1 4 



|0) 



(0| 



6T 
j7(e) 



1 2 



|0) 



IV. CLASSICAL MODEL FOR THE STATISTICS OF THE MICROWAVE 

BACKGROUND ANISOTROPIES 

A distribution of the microwave background temperature over the sky is a real function of 
the angular coordinates. Assuming that 6T/T is a sufficiently smooth function on a sphere, 
one can expand it over the set of orthonormal complex spherical harmonics Yim{0, 0) [15]: 

JTrj-i 00 I 

^(e) = E E hmYimie) + a^^rrje)] . (30) 

We want to formulate a statistical hypothesis about the coefficients aim, so it is better to 
write them first in terms of real (r) and imaginary (i) components: 

00 I 

2E E WlmYLi^) - O^lmYLi^)] ■ (31) 
l=Q m=-l 

Our statistical hypothesis is as follows: (i) all members of the set of random variables 
{a[^, a]^} are statistically independent, (ii) each individual variable is normally distributed 
and has a zero mean, (iii) all variables with the same index I have the same standard 
deviation it;. All said is expressed by the probability density function (p.d.f.) for individual 
variables: 



O'lm — 



ST 
T ^ 
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and by the p.d.f. for the entire set of variables, which is simply a product of all p.d.f.'s for 
all individual variables: 

/ = ni,m/(aL)/(a^„) . (33) 

Having postulated the p.d.f.'s, we can now compute the expectation values of certain 
functions of the random variables. Below, the angular brackets will denote the expectation 
values calculated with the help of the p.d.f. (33), unless other definition is stated. 

Obviously, all linear functions have a zero mean: 

(34) 

. (35) 

All triple products have zero means. Among quartic combinations, only those can survive 
which have four indices (r), or four indices (i), or two indices (r) and two indices (i). Two 
representative expressions are: 

^'^hmi^hm2^hm3^km4^ ^ '^h(^f3^hl2^mim2^l3l4.^m3m4 
+ Crh<^l2^lil3^mim3Sl2hSm2rrn+ (^l^(^l2^hlJmim4Sl2l^Sm2m3 , (36) 
<(^hmi(^hm2^hm3^km4> ^ (^h(^l3^hl2^rnim2^l3lJrn3m4 ■ (37) 

Other quartic combinations can be obtained by the replacement (r) (i) in Eqs. (36), (37) 
(or by permutation of pairs (Im) in case of Eq. (37)). The higher-order correlations can be 
derived in a similar way, but we will not need them. 

In our further calculations related to the random variables ^(e) and ^(ei)^(e2) it is 
easier to deal with the complex coefficients aim, so we will first translate the above relation- 
ships to them. By using the available information one can derive 

<04rn> =0, < 0^hmi(il^m2^ = '^Cff 5l^l25mim2 i 
< Ct/imi<2/2m2<^/3m3'^Um4 > ~ ^^h^'l2i^hl3^m'i_m35l2li^m2m4, + ^hU^m'i_mi5l2l3^m3m3) ■ (38) 

The mean values of the complex conjugated quantities are given by the same formulas (38). 
Other nonvanishing quartic combinations can be obtained from the one in Eq. (38) by the 
permutation of pairs {Im). 

Now, even before deriving the p.d.f.'s for the random variables ^(e) and ^(ei)^(e2), 
we can find some expectation values. It is clear from the definition (30) and Eq. (38) that 

(^(e)^=0 . (39) 

When calculating the angular correlation function one should remember that 



For quadratic combinations we have 

T" T' 2 C C" I Z 

<^hmi^l2m2 > ~ 0-lj^\l2^mim2 i < (^Iimi^l2m2 > 
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' 2/ + 1 

E ll^(ei)r4(e2) = ^PKcos5) 

m=—l 



(40) 



(note the origin of the factor 2/+ 1 which will accompany us often). By taking the product of 
two expressions (30) and using Eqs. (38), (40) one can find the angular correlation function 



^(ei)^(e2)^ = -E^f(2^ + l)PKcos5) 



(41) 



If the separation angle 5 is zero, we obtain 

(f(e)f(e); 



1 oo 
^ 1=0 



(42) 



[One may notice an incidental fact that the mean value of the random variable af defined as 



is <a'i>-- 



2(2/ + 1)(T; , that is the same expression which enters Eq. (41). 



This may suggest an interpretation of the quantity <af> — <af>'^ as the variance of the 
multipole moments. One should be carefuU with this interpretation, see Sec. V.] 

We can also find the 4th order expectation values. The product of 4 expressions (30) in 
conjunction with Eq. (38) gives 



'5T ST ST ST \ /ST ST ' 
^(ei)^(e2)^(ei) — (62) ) - ( — (ei) — (eg) 

/ST^ ST^ X IST ^ ST^ y 
= (— (ei) — (62)) +(— (e) — (e)\ . 



(43) 



If 5 = 0, it follows from Eq. (43) that 

1 4 



ST 



e) 



ST 



2v 2 



(44) 



Up to difference in the meaning of the angular brackets, the formulas (39), (43), (44) 
reproduce the analogous results of the previous Section. Moreover, from comparison of 
Eqs. (23), (24) with Eqs. (41), (42) we can relate the quantities Ki, derivable from a given 
cosmological model plus perturbations, with the abstract quantities u/. 

We can now engage in our major enterprise — the construction of the p.d.f. for the 
random variable v = ^(61)^(62). We will start from the p.d.f. for the random variable 
z = Y'i^)- When it is necessary to distinguish directions ei and 62, we will use the notations 
Zi and Z2. 

The variable ^ is a function of the variables {a[^, a]^} whose p.d.f. 's are known, Eqs. (31), 
(32). There exist regular methods (see, for example, an excellent book [16]) allowing to derive 
rigorously the p.d.f. of a function. However, in our case that the function is linear and all 
p.d.f 's are normal, we can partially rely on a guesswork. Combining formulas and guessing 
we can write 



1 



27r (72 



(45) 
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where 

1 oo 

^' = -E^'(2^ + l) ■ (46) 

^1=0 

The p.d.f. (45) certainly leads to Eqs. (39), (42), (44). Moreover, it allows us to say that 
the probability to find z outside of la^ interval is approximately 32%: 

P{\z\ > a^) ^ 0.32 . 

We now introduce two variables, Zi and Z2, and ask about the p.d.f. in the 2-dimensional 
space {zi, Z2). Again, partially relying on a guesswork, we find that 

/(-I, = exp + 4 - 2pzrZ2]^ (47) 

where 

P^'=-E(2^ + lW^'(cos5), |p|<l . (48) 

^1=0 

(See Eq. (5.11.1) in Ref. [16]). First, we can check that the marginal distributions are correct. 
For f{zi), one obtains 



/oo _ 
f{zi,Z2)dz2^ r— e 2^ 
-oo Y ZTT Gz 



and one obtains a similar expression for f{z2)- Second, one can check that 

{zD = al , l^zl) = al , {Z1Z2) = pal 

where the angular brackets mean the integration with the p.d.f. (47). These equalities are 
Eqs. (42), (41) which we must have obtained. 

Finally, wc shall derive the p.d.f. for the variable v = ZiZ2- We will do this in some detail 
following the prescriptions of [16]. 

Let us introduce the two new variables {zi, v) instead of {zi, Z2) according to the trans- 
formation 

V 

Zi = Zi , Z2 = — 

Zl 

The Jacobian of this transformation is J — 1/zi. The p.d.f. f{v) is the result of the 
following integration: 



/(^) = 7^ o \ -> / r-l exp 



27ra2^/r^y-oo l^il ^[ 2(T2(l-p2 



2 V ^ 

Zi + ^ - 2pv 
zi 



dzi . 



The integral over zi can be taken with the help of 3.471.9 from [17]. The resulting p.d.f. 
can be written in the form: 
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f{v) = 



pv 

-^e^p') K. ( -271^ ) , for ^; > 



i ^^^^'^^o (^) , for.<0 



where Kq is the modified Bessel function of its argument [18]. 

The function f{v) is quite complicated and the distribution is obviously not normal. The 
function f{v) goes to zero for v — > ±(X) and diverges logarithmically at the point v — 
Even a verification of the normalization condition 

/oo 
f{v)dv = 1 (50) 
-oo 

is not trivial. However, with the help of 6.621.3, 9.131.1, 9.121.7, 1.624.9 and 1.623.2 
from [17] one can prove the validity of Eq. (50). 

The mean value and the standard deviation of the variable v are known, see Eqs. (41), 
(43): 



{v) = pal , <7„ = [{v') - [vf] = al4^l . (51) 

We already knew that the standard deviation is big. We now see again that a^ = a^ at the 
separation angles at which the angular correlation function vanishes, p — 0, and a^ — V^a^ 
at zero separation angle, p — 1. For other separation angles, a^ lies between these two 
numbers. 

Now that we know the p.d.f., we can assign probabilities to the different ranges of the 
variable v. For instance, we can calculate the probability that the measured v will be found, 
say, outside of the A(T„ interval surrounding the mean value of v, where A is an arbitrary 
fixed number. The probability of our interest is 



P{\v-{v)\ > \a,) = / f{v)dv+ ^f{v)dv . (52) 

To get a qualitative estimate of the associated theoretical uncertainties for the observable 
V, we will ask a slightly different question. What should the number A be in order to have 
the 0.32 chance of finding v outside the Xa^ interval and, hence, the 0.68 chance to find it 
inside the interval? 

To evaluate the size of the disaster, we will start from the case p = 0. In this case, the 
p.d.f. (49) is symmetric with respect to the origin v — (this is why (v) is zero in this case) 
and 

P{\v\ > Xal) = - / Ko{x)dx . (53) 

We want this number to be approximately equal to 0.32. Judging from the Fig. 9.7 in 
Ref. [18], a half of the area under the Kq{x) function is accumulated when integrating from 
approximately a; = 1/2 and up to infinity. This means that A should approximately be equal 
to 1/2 . 

If p 7^ the evaluation of P is more complicated. For p ^ 0, the function (49) is not 
symmetric with respect to the origin v — 0. It has larger values at positive v's if p > (this 
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is why (y) > in this case) and it has larger values at negative v^s if p < (this is why 
(v) < in this case). The graph of the function e^Ko^x) plotted on Fig. 9.8 in Ref. [18] 
is helpful. A qualitative analysis shows again that A is approximately equal to 1/2. (More 
accurate estimates can of course be reached by numerical methods.) 

At any rate, the |(T^ interval gives approximately the same probability estimates as if 
the distribution (49) were normal. 

V. ON THE "COSMIC VARIANCE" 

The set of random variables {a[„,, ci^m} defined by Eqs. (32), (33) lives its own independent 
life regardless of whether or not the variables are considered random coefficients in the 
expansion of some function over spherical harmonics. Being such, it allows introduction of 
new functions and calculation of their expectation values. One interesting variable is defined 
by the equation 

I I I 

m=—l m=—l m=—l 

By using Eq. (38) one can calculate the expectation value of a^: 

{af)^{2l + l)2af . (55) 

The factor 2(2/ + 1) reflects the number of independent "degrees of freedom" associated with 
the index /. One can also introduce the variable af and calculate its expectation value: 

(af) = {21 + m + l)8af = {ajr^^±^ . (56) 

The difference (af) — {afy is, by definition, the variance of the variable af. Prom Eqs. (56), 
(55) one finds 

(4) - {ah' = ^{ah' . (57) 

This formula, as it stands, expresses a well-known fact: the variance of the random 
variable defined as the sum of squares of n independent random variables (degrees of 
freedom) with the same normal density, is | times smaller than the square of the mean 
value of [16]. There is nothing "cosmological" or "infiationary" in this fact. Knowing 
the p.d.f.'s for the set {a[^, aj^} one can calculate the higher-order correlation functions 
for the variable af and its distribution function [23-29]. In the recent literature, formula 
(57) became known as the "cosmic variance" . Formula (57) and a possibility (or lack of) to 
extract complete information about a stochastic process from its single realization are, in 
general, different issues. For ergodic processes, the existence of a definitely true relationship 
(57) prevents in no way the extraction of complete information about the process from a 
single reahzation [30]. 

It is important to realize that it is the mean value of the random variable af, not af itself, 
that enters the expected angular correlation function in front of the Legendre polynomials 
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(54) 



and which is often called the multipole moment. The af is a random variable and its variance 
has meaning, the (af) is a number and its variance has no meaning. Specifically, one can 
notice that the angular correlation function (41) can be written in the form 



ST ST \ 1 °° 



T 



27r 



(58) 



On this ground, there may be a temptation to write the random variable ^'{^1)^(^2) in 
the form 



ST. .ST. . 1 ^ 

^(ei) — (62) = — ^afP^(cos(5) 



(59) 



and to interpret Eq. (57) as the variance for the multipole moments of the correlation 
function. One should resist to this temptation. 

Let us show that the definition (59) is incorrect despite the fact that it gives correct 
expectation value (58). It follows from the definition (59) that 



ST ST ST ST 
— (ei) — (62) — (ei) — (62) 

-i 00 1 00 

= r^E«f[^Kcos5)]2 + — j2 a]Pi{cosS)alP,{cosS). 

^"^ 1=0 ^'^ l,l'=0,l^l' 



(60) 



Using (56), (58) and remembering that af and af are statistically independent for I ^ I', 
one can find the expectation value of the quantity (60): 



'ST ST ST ST ; 
— (ei) — (62) — (ei) — (62) 

1 °° 1 
E(af)[PKcos5)]^ + 



47r2 



1=0 



47r2 



£ {af){af)Pi(cosS)Pu{cosS) 



l,l'=0,l^l' 

00 

J2{af)P,icosS) 

Ll=0 



-1 00 

£-,u^fnPi{coss)r 
1=0 



'ST ST 
^(^1)^(62) 



00 

/=0 



1 



47r2 2i + 1 



{afy[Pi{cosS)Y. 



(61) 



It follows from (61) that the variance of the variable ^(61)^(62) would read (if (59) were 



correct) : 



-1 00 

-Y,{2l + l)anPi{cosS)f 



TT 



(62) 



1=0 



This expression should be compared with the correct variance following from Eq. (43): 



1 



1 2 



J2{2l + l)afPi{cosS) 

Ll=0 



1 



1 2 



E(2^ + ^)^f 

.1=0 



(63) 



Formulas (62), (63) disagree even for S — 0, and even in their first, / = term. This shows 
that the ad hoc definition (59) is incorrect. The correct definition of the random variable 
^(61)^(62) is the one following from the definition (30) and which we have used in this 
paper. 
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VI. CONCLUSIONS 



A particular cosmological model plus perturbations gives unambiguous predictions with 
regard to the expectation values of the measurable quantities. Differing models give different 
predictions. We want to distinguish them observationally and to learn about physics of the 
very early Universe. However, the quantum-mechanical origin of the cosmological pertiir- 
bations is reflected in the theoretical statistical uncertainties surrounding the expectation 
values. One important measurable quantity is the angular correlation function of the mi- 
crowave background anisotropies. Its mean value at the zero separation angle was denoted 

in this paper. It was shown that the standard deviation for the correlation function is 
very big. The 68% confidence level corresponds, approximately, to |(T^ at the separation 

angles where the correlation function vanishes, to -^cr^ at the zero separation angle, and to 
intermediate numbers for other separation angles. 

The angular correlation function has actually been measured. It is presented at Fig. 3 
in the paper [19]. The authors surround the measured points by a narrow shaded region 
which they address as follows: "The shaded region is the 68% confidence region .... including 
cosmic variance and instrument noise" . It is not quite clear what the authors of Ref. [19] 
(see also Ref. [20]) mean by "cosmic variance", but if they mean the theoretical statistical 
uncertainties for the correlation function variable v, these uncertainties are significantly 
larger than what is plotted. According to the calculations presented above, the half- width of 
the shaded region should be approximately 600 {jJ-KY near the points where the correlation 
function vanishes and approximately 840 {/J'K)^ near the point marking the zero separation 
angle. 

The conclusion is a bit disappointing. Apparently, God is telling us something important 
about the very early Universe by exhibiting the microwave background anisotropies, but the 
channel of information is so noisy that it will be hard to understand the message. 
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APPENDIX. ON THE "STANDARD" INFLATIONARY FORMULA FOR 

DENSITY PERTURBATIONS 

In the body of this paper we have studied the statistical properties of cosmological pertur- 
bations of quantum-mechanical origin and related anisotropics in the microwave background. 
These properties arc essentially universal, they arc equally well true for perturbations of any 
nature - density perturbations, rotational perturbations, or gravitational waves. However, 
in practice, it is very important to know which sort of perturbations we are actually deal- 
ing with, assuming that the observed large-angular-scale anisotropy is indeed caused by 
perturbations of this origin. In other words, what does theory say about the comparative 
values of the amplitudes, if we agree on the generating mechanism, equations, and a class 
of cosmological models of the very early Universe, say, governed by the scale factors (15)? 
In addition to gravitational waves, density perturbations can also be generated by the same 
mechanism, if one makes favorable assumptions about the dominant matter in the very early 
Universe (scalar field) and its coupling to gravity (minimal, the same as for gravitational 
waves). According to calculations of Rcf. [3], density perurbations and gravitational waves 
will have amplitudes of the same order of magnitude, wheras according to the "standard" 
inflationary formula the amplitude of density perturbations will be many orders of magni- 
tude larger than the amplitude of gravitational waves, if the expansion rate of the very early 
Universe was sufficiently close to the archetype inflationary model - the de Sitter expansion. 

It is necessarry to say that the quantum-mechanical generating mechanism has become 
very popular in the context of the inflationary hypothesis. Inflationary literature often speaks 
about cosmological perturbations being generated "from quantum fluctuations" . However, 
this literature associates the explanation of the phenomenon with such things as ambiguity 
in the choice of time in the de Sitter universe, horizon temperature, tremendous inflation of 
scales, and so on. The basic concepts are adjusted accordingly. Instead of amplification, with 
the emphasis on a nonvanishing parametric coupling, increase of amplitude at the expense 
of energy of the pump fleld, quantum-mechanical generation of waves (particles) in strictly 
correlated pairs, etc., inflationary literature speaks about magnification, with the emphasis 
on "stretching the waves" and "crossing the horizons" . Apparently for these reasons, infla- 
tionists did not get puzzled with their "standard" formula for density perturbations which 
states that one can produce arbitrarily large amount of density perturbations by practically 
doing nothing. 

The "standard" formula relates the amplitude of density perturbations today with the 
values of the scalar fleld during inflation. Let us consider, for deflniteness, perturbations of 
the matter density Sp/p with today's wavelengths of the order of today's Hubble radius Ih- 
The "standard" formula says that 



6p 
P 



(64) 



where the right hand side of this formula is supposed to be evaluated at the time ti when 
the wavelengths of our interest were "crossing the horizon" during inflation. Let us agree 
with the so-called "slow-roll" approximation and assume that the Hubble parameter H was 
almost constant during that epoch, \H\ <^ H^. Let us take the numerical value of H during 
that epoch at the level, say, 20 orders of magnitude smaller than the Planck value of H. 



22 



For quantum-mechanically generated gravitational waves, this would result in today's am- 
plitude h ^ 10~^° and the induced anisotropics of the microwave background 5T /T ~ lO^^*' 
which are much much lower than the level currently discussed in the experiment. However, 
for density perturbations, according to the "standard" inflationary formula, the situation is 
totally different. Without changing anything in the curvature of the space-time responsible 
for the generating process (that is, leaving H almost constant and at the same numerical 
level 20 orders of magnitude smaller than the Planck value), but simply sending (f>{ti) to zero 
(which corresponds, due to the Einstein equations, to sending H to zero, i.e., making the 
"slow-roll" approximation better and better, making the expansion law closer and closer to 
the de Sitter expansion and making the shape of the generated spectrum closer and closer 
to the scale-invariant form) one produces arbitrarily large {Sp/p)\H- Inflationists love to 
stress that the de Sitter gravitational pump field generates perturbations with the Harrison- 
Zeldovich (scale-invariant, flat) spectrum. What they do not stress is that, according to the 
"standard" formula for density perturbations, the amplitudes of the scale- invariant spectrum 
are infinite, and the amplitudes of the almost scale-invariant spectrum are almost infinite. 
Instead of blaming their own formula, inflationists blame the scalar field potentials. This 
formula is the reason for rejecting certain scalar field potentials on the grounds that they 
generate "too much" of density perturbations, for claims that the contribution of gravita- 
tional waves to ST/ T is "negligibly small" in the limit of the de Sitter expansion, and even 
for claims about copious production of black holes during inflation. 

Recently, the "standard" formula has been claimed to be confirmed [31] and recon- 
firmed [32]. In Ref. [31], this formula has been formulated, essentially, as the following 
"standard result": "... we see that the scalar perturbations can be very strongly amplified" 
[the increase of numerical value from (almost) zero to (almost) infinity] "in the course of the 
transition" [the instantaneous change of the cosmological scale factor from one power-law 
behavior to another power-law behavior]. The authors of the paper [31] assure the trusting 
reader: "We think that there is nothing strange about this ...". 

Here, we will try to understand the origin and mathematical justification for the "stan- 
dard" infiationary formula. 

The early papers which are usually quoted in this connection are the papers [33-35] . We 
will start from the paper of Hawking [33] which seems to be clearer than others in expressing 
the basic idea and intentions. The papers [33-35] are similar in many respects. 

Hawking considers a scalar field (6 running slowly down an effective scalar field potential. 
He discusses the inhomogeneous fluctuations 0i(t, x) in the fleld (p = (j)o{t) + 0i(t, x) which 
mean that on a surface of constant time there will be some regions where the (p field has 
run further down the hill than in other regions. He introduces a new time coordinate 
t = t + 6t{t, x) in such a way that the variations of the field are removed and the surfaces of 
constant time are surfaces of constant (p. Since the scalar field transforms as 
00 + 01 ~^ 00 + 01 — 00^^; the required condition is achieved by the time coordinate shift 6t = 
0i/0o. Then Hawking says that the change of time coordinate will introduce inhomogeneous 
fluctuations in the rate of expansion H. He and other authors take (apparently, on the 
grounds of dimensionality only) SH ~ H^dt. From here they come, implicitly or explicitly, 
to the dimensionless amplitude of density perturbations 



5p SH 



— ~ HSt 



(65) 



P H 
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Some authors write explicitly (f)i ^ H and 5p/ p ~ if^/0o- 

The analysis has been done at the inflationary stage. To obtain the today's amplitude 
of density perturbations in wavelengths, say, of the order of the today's Hubble radius, it is 
recommended to calculate the right hand side of Eq. (65) at the moments of time when the 
scales of our interest were "crossing horizon" during inflationary epoch. In one or another 
version this formula appears in the most of inflationary literature and because of numerous 
repetitions it has grown to the "standard" one. According to this formula, the amplitude of 
density perturbations becomes larger if one takes the <^o smaller. 

The authors of [33-35] work with a speciflc scalar field potential, so the numerical value of 
00 and the numerical value of 5p/ p following from Eq. (65) turn out to be dependent on the 
self-coupling constant in the potential. These authors are concerned about the unacceptably 
large amplitude of density perturbations that they have produced. But this is not a concern 
about the fact that the Einstein equations play no role in this argumentation, it is a speciflc 
detail in the scalar fleld potential that the authors of [33-35] do not like. 

Now let us show the shortcomings of the argumentation in [33-35]. Let us consider a 
scalar fleld (f) with arbitrary potential. Write the fleld as = 0o(O + 'Pi{'t)Q where Q is the 
n-th spatial harmonic, Q'^^ + r?Q — 0. Write the perturbed metric in the form 

ds^ = -dt^ + a^{t)\{\ + h{t)Q)5ij + hi{t)n-^ Q,i,j]dx' dx^ . 

The de Sitter solution corresponds to (f)Q = 0, a{t) ~ c^*, and H{t) = a/a = const. It follows 
from the Einstein equations that the (linear) contribution of the scalar fleld perturbations 
to the total energy density e — to + can be written as 

e</> = 00 {^i - 0i[ln(a^0o)]'} Q ■ 

The contribution e^, as well as other components of the perturbed energy-momentum tensor, 
vanish in the de Sitter limit 0o ~^ 0. 

Thus, the flrst conclusion we have to make is that in the de Sitter limit there is no linear 
density perturbations at all. The scalar fleld perturbations are uncoupled from gravity, they 
are not accompanied by linear perturbations of the energy-momentum tensor and they are 
not accompanied by linear perturbations of the gravitational fleld. The general solution 
to the perturbed Einstein equations is a set of purely coordinate solutions which can be 
produced or totally removed by appropriate coordinate transformations. The scalar fleld 
perturbations reduce to a test fleld whose only role is to identify events in the spacetime. 
One can still ask about a coordinate system such that the surfaces of constant time r, 
T = 0i(t,x) are surfaces of constant cf). But the perturbation of the expansion rate of this 
new coordinate system will have nothing to do with perturbations in the energy density. 
Despite the presence of the test scalar fleld, every space-like hypersurface is a surface of 
constant energy density. 

Now let us assume that (po is not zero. Transformation of time i = t + x{t)Q generates 
a Lie transformation of the scalar fleld: 

Mt) + Mt)Q ^ Mt) + IMt) - Mt)x{t)]Q ■ 

If one wants the transformed fleld to be homogeneous one takes x(^) = <f>i{'t) / 4>o{t) ■ The 
same transformation of time generates Lie transformations of the metric. The transformed 
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Qoo component is cjoo = — 1 + SxQ, the transformed Qik components are described by 
h — h — 2{a/a)x- There appear also the components but they will not participate in our 
linear analysis. The expansion rate of a given frame of reference is determined by the trace 
of the deformation tensor [36] : 

In the linear approximation and before the transformation, 

D ^ 3H + hsh - hi)Q . 

After the transformation, D = D — 3HxQ- 

So, the introduced inhomogeneous fluctuation in the rate of expansion is 5H — —HxQ = 
H5t, not 5H — H'^St assumed in Refs. [33-35] . 

The Einstein equation for energy density D^/3 = Ke is satisfied before and after the 
transformation, since the variation of H is balanced by the variation of the energy density. 
The transformed energy density is 

e = €o + €^- ioxQ = eo + 00 ( T" ) Q ■ 

\(PoJ 

Thus, if one makes the 0o smaller, the energy density perturbation decreases according 
to the Einstein equations, and it increases according to the conjectures of Refs. [33-35]. The 
use of the correct expression SH ~ H6t in Eq. (65) makes the Sp/p decreasing when the 0o 

is decreasing. 

The situation becomes even more disturbing if one recalls that the formula (64) has been 
seemingly confirmed and derived rigorously as a result of more detailed studies. People did 
really write the perturbed Einstein equations. Moreover, it was done in the framework of the 
so-called gauge-invariant formalism, the whole purpose of which is to eliminate coordinate 
solutions and to work exclusively with something "physical" . The basic mathematical tool 
in these studies is the gauge-invariant potentials $ and "if constructed from the components 
of the perturbed metric. 

As we have seen above, density perturbations in the scalar field matter vanish when 0o 
goes to zero, and there is no density perturbations at all at the de Sitter stage. This is true 
irrespective of the wavelength of the perturbation. It can be shorter or much longer than 
the Hubble radius, that is, it can be "inside" or "outside" the Hubble radius. In particular, 
this is true of the perturbation whose wavelength is such that it will grow by today to the 
scale of today's Hubble radius. The gauge-invariant potentials $ and ^ are strictly zero at 
the de Sitter stage. Why does then the today's amplitude of the perturbation go to infinity, 
according to the "standard" formula, in the limit 0o — at the "first horizon crossing"? 
Because, the infiationary literature explains, the amplitude will be almost infinitly enhanced 
in "the course of transition" of the background equation of state from the quasi-de Sitter 
one p ~ — e to the radiation-dominated p = |e or matter-dominated p — one. The 
favorite concept in this argumentation is the "constancy of (" . The often quoted papers are 
Ref. [37], which uses notations and equations of [38], and Ref. [39] which summarizes the 
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previous work and gives a clearer exposition. We will follow equations and notations of the 
paper [39]. The advantage of this paper is that it contains enough mathematical details to 
make it possible to follow the spirit and the letter of calculations. 

Mukhanov et al. [39] work with the gauge-invariant potentials $ and ^. According to 
one of the perturbed Einstein equations, $ = ^f. In terms of the the basic equation of [39] 
at the scalar field stage is 

+ 2^^^$' _ + 20[, f 5 V * = (66) 

(«/0o) \n) 

where ' = d/dr), dt = adr], 7i = a' /a. Equation (66) is exactly the same equation as the 
basic equation (2.23) of Ref. [37]. At the perfect fluid stage, the basic equation of [39] for 
"adiabatic" density perturbations is 

$" + 37i(l + c^)$' - c^V^$ + [2n' + (1 + 3cl)n'^]^ = (67) 

where = Sp/Se = po/^o- Equations (66) and (67) were derived from the original perturbed 
Einstein equations with the help of manipulations aimed at expressing the equations in terms 
of the gauge-invariant potentials. 

The authors of [39] introduce a new quantity C. defined as 

where w = Pq/^q. This quantity is simply a new letter. As soon as the function $ is known, 
the function ( can be calculated from the definition (68). Using the definition of Eq. (66) 
can be written in the form 



and Eq. (67) in the form 



^CH{l + w)^\v'^ , (69) 



^CH{l + w) = \clV'^ . (70) 



Mukhanov et al. [39] consider perturbations with wavelengths "far outside the Hubble 
radius for which can be neglected". Neglecting the right hand sides of Eqs. (69), (70) 
they arrive, in this approximation, at the equation C ~ and the "conservation law" 

C const . (71) 

It is important to note that the derivation of this conservation law did not require any 
knowledge about the initial data and solutions for $. The constant in the right hand side 
of Eq. (71) emerges as a universal number, irrespective of any particular solution for $. 
Although the "constancy of is a favorite notion in the inflationary literature, it appears 
that inflationists have never asked what the origin and numerical value of this constant is. 
If this constant is a universal number, is it equal to a billion, or one, or zero? We will later 
show that this particular constant must be equal to zero. However, without specifying the 
value of this constant, it is often assumeed that it is not zero. 
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The next step in this argumentation proceeds as follows. Imagine that the background 
equation of state changes during a short interval of time from the initial (i) p ^ —e to the 
final (/) p = |e or p = 0, so that 1 + Wi ^ 1 while 1 + ~ 1. The authors of [39] consider 
the evolution of a long- wavelength perturbation from the initial Hubble radius crossing (tj) 
to the final Hubble radius crossing (tf). They make additional assumptions, they assume 
that $ vanishes "at very early and very late times". In the definition (68), they drop the 
term with $ and simplify the quantity 



3(l + 'y;(t/)) 

Then they refer to the constancy of C, C(^i) = C(i/)) ^-iid arrive at the formula 

l + w(ti) 5 + 3w(t/) l + w(ti) ^ ' 

(which, in fact, is in a conflict with the constatncy of C, as we will show later). 

According to this formula, and since 1 + w(t^ can be arbitrarily close to zero, the 
$(t/) can be made arbitrarily large for any nonvanishing $(tj). Moreover, since the time of 
transition from one equation of state to another can be arbitrarily short, the tremendous 
increase of numerical value of the potential $ is supposed to happen almost instantaneously. 
(The author of [32], who reconfirmes the "standard" results, has even plotted a graph for 
this jump.) The authors of Ref. [39] emphasize that their result (72) (the actually published 
[39] formula (6.67) contains a misprint: the position of symbols w{ti) and wi^j) should 
be interchanged) is in a full agreement with previous studies and Eq. (64). Formula (72) 
suggests an arbitrarily large production of density perturbations for no other reason but 
simply because the 1 + w(tj) was very close to zero. There is something strange with this 
formula. [I realize well that what I qualify here as strange is certainly considered by others 
as perfectly alright. Otherwise somebody would raise a voice of protest against the ease 
with which inflationists generate tremendous amounts of various substances (some of the 
authors are even claiming that they can "overdose" our Universe). However, judging from 
the literature, it is not only that there are no voices of protest, but there is rather an element 
of competition as for who was the first to proclaim the "standard" infiationary results. For 
instance, the authors of [40] address the inflationary claims about density perturbations as 
"flrst quantitatively calculated in [33-35] [and which] have been successfully quantitatively 
conflrmed by the COBE discovery".] 

Let us now try to sort out what we are dealing with. 

The flrst point to realize is that the basic equations (66), (67) are, strictly speaking, 
incorrect. They are incorrect in the following sense: they are constructed from a correct 
equation and the flrst time derivative of the correct equation. Transforming the original 
Einstein equations, the authors of [38,39] have effectively raised the order of differential 
equations. If the correct equation is satisfled, the constructed equation is satisfled too, 
but not vise versa. To get the feeling of the danger involved, consider a correct equation 
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X — and a constructed equation x^Ax = where A is arbitrary constant. Solutions to the 
correct equation do satisfy the constructed equation, but the latter one admits exponentially 
growing solutions which are not allowed by the original equation. 

Let us start our analysis from Eq. (66). Since the potential $ is an eigenfunction of the 
Laplace operator, = — we will replace with — We introduce also a new 
function of the scale factor a{rj): 

7 = -H/H^ = 1 + {a/a'y 

and use the background Einstein equations in order to express the coefficients of Eq. (66) in 
terms of the scale factor a{r]) and its derivatives. Now, introduce a new variable according 
to the definition 



1 a' 



-7 



277,2 a y 



(73) 



So far, the variable fi is simply a new variable replacing $, but the importance of fi is in 
that the original perturbed Einstein equations require this variable to satisfy the equation 



2 («V^)" 

n — 







(74) 



(For those interested in gauge-invariant potentials, I may remark that is a genuine gauge- 
invariant variable in the sense of [38,39].) With the help of Eq. (73), Eq. (66) identically 
transforms to 



0. 



(75) 



If Eq. (74) is satisfied, Eq. (75) is satisfied too, but not vice versa. Equation (75) is totally 
equivalent to 



0^7 



a^7 



/i 



+ n' 



«V7 



X 



(76) 



where X is arbitrary constant. Thus, Eq. (66) involves an arbitrary constant X which, in 
fact, must be zero according to Eq. (74). 

Let us now turn to the perfect fluid equation (67). The situation here is similar to the 
scalar field case. Introduce a new variable v according to the definition 



1 a' 7 
2n2 a ci 



(77) 



The importance of u is in that the original perturbed Einstein equations require this variable 
to satisfy the equation 



2 2 

n cz 



= 



(78) 



28 



(The function z/ is a genuine gauge-invariant variable.) With the help of Eq. (77), Eq. (67) 
identically transforms to 




= 



(79) 



This equation is totally equivalent to 



[ \a 




V 




a 



Y 



(80) 



where Y is arbitrary constant. Thus, Eq. (80) does also involve an arbitrary constant Y 
which, in fact, must be zero according to Eq. (78). The role of these constants X and Y in 
the constancy of ( argument we will discuss shortly. 

Equations (74), (78) is all we need to solve, in order to find all the perturbed components 
of the metric tensor and energy-momentum tensor. There is no other way to find the 
evolution of density perturbations except of doing the hard work of solving these equations, 
imposing appropriate initial conditions, joing the solutions, etc. (This is what is being done, 
convincingly or not, in Ref. [3]. It is only great physicists can afford having one or several the 
greatest mistakes of their life, wc can not afford any.) Among other things, these solutions 
will tell you what is the value and time dependence of the functions C introduced by the 
definitions (68), (73), (77). The correct equations do not show anj^thing like enormously 
large amplitude of today's density perturbations in the limit of the de Sitter expansion. 
However, we need to return to the concept of "constancy of (" which pretends to answer 
important physical questions without solving any equations at all. 

Combine the definitions (68), (73) to show that ( at the scalar field stage is 



The term with the Laplacian neglected in Eqs. (66), (67) is the term with in equations 
(76) and (80). Let us drop this term as the authors of [39] do. Then, in this approximation, 
Eq. (76) gives ( ^ X/2n'^ — const and Eq. (80) gives ( ^ Y/2n'^ — const. These relationships 
show that C can be a universal constant, independent of any particular solution, only if it 
is supported by the constant X or, correspondigly, by the constant Y. But, as we have 
shown above, these constants must be equal to zero, if one is willing to work with correct 
equations. The "constancy of argument fails at the very beginning, regardless of validity 
or not of the additional assumptions that have been made on the route to Eq. (72). The 
conservation law ({ti) — C(^/) degenerates to an empty statement = 0, and nothing can 
be derived from it. 




Combine the definitions (68), (77) to show that ( at the perfect fluid stage is 
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The idea of a free lunch based on the "constancy of (" seems to be so attractive that 
the following question is often being asked. Suppose that not for all solutions for $, but for 
some of them, suppose that not in the leading order n~^, but in the next order, suppose, 
nevertheless, that ( calculated from this specially chosen solution for $ is approximately a 
nonzero constant (which is indeed possible). Why cannot we return to Eq. (68) and repeat 
all the arguments that have seemingly led us from the definition (68) to the result (72)? 
For instance, the authors of [31] have even considered, along this line, a concrete model 
consisting of two consecutive power-law scale factors. They speak about an initially small 
potential $ being "distributed" after the transition point "into two modes, both very large 
in amplitude, one which decays, the other yielding" the "standard" result (72). And this is 
how, they argue, the "scalar perturbations can be very strongly amplified in the course of 
the transition". 

So, we also need to consider this model and explore what is contained in the definition 
(68) treated as an equation for 

First of all, use the background equation 



2 H 



and write C in a more convenient form: 



aH \H 

Integrate this equation to produce the general solution 



a 



C + 



10 



dt 



(81) 



(82) 



where C is arbitrary integration constant. So far, ( can be an arbitrary function. Now 
assume that C is a constant, C = Co = const, and write 



* = Co 



H ft 
a Jti 



adt 



a 



(83) 



[The constants Co and C can be related to the coefficients in front of two linearly independent 
solutions to Eq. (74). In the long wavelength approximation. 



1 



n 



a^7 



a^'^dr] ) dr] 



a^7 



Using the definition (73) one can show that Co = ^ = 

As the second step, let us consider, together with Deruelle and Mukhanov [31], an expand- 
ing model which decribes a transition from one power-law scale factor to another power-law 
scale factor. Let us write 



ttf = a2{t — 



\P2 



(84) 



The scale factor and its first time derivative are continuous functions at the transition point 
t — ti. This requires 
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P2 



To be closer to the inflationary results, one can keep in mind that pi ^ 1, while p2 is 1/2 
or 2/3. 

Now we will follow the time evolution of $ specifically in this model. Let tj be the time 
when our perturbation first crosses the Hubble radius: a(ti)H{ti) = n. Let be the time 
when the perurbation returns back inside the Hubble radius: a{tf)H{tf) = n. The functions 
a{t) and H{t) are continuous all the way from ti to tf, and so is the function (83). The 
initial value of the potential is 

m)^Co + C^ , (85) 



the final value is 



a{ti) 



^tf)^CoI + C^ (86) 



a{tf) 



where 



a(tf) Jti 



Combining (85) and (86) we find 



adt 



H{tf) HjU)^ 
_a{tf) a{ti) 



(87) 



The integral / can be easily calculated for the scale factor (84). Since a{tf) ^ a(ij), pi ^ P2, 
and a{ti)H{ti) ^ n, we have approximately 



P2 + 1 

Formula (87) reduces to 



'P2 + 1 a{ti) P2 + 1 

The constant C could be set to zero from the very beginning. In any case, the term 
C {H (ti) / a{ti)) is of the same order or smaller than $(ti). By setting C = 0, we arrive 
at 

P2 + ^ 

There is nothing like tremendous jumps of $ at the transition point. The "constancy of (" 
effectively translated into constancy of [The fact that $ is approximately constant for 
wavelengths longer than the Hubble radius follows also from the exact solution to Eq. (74) 
which can be found in case of power-law scale factors.] 
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The best what we can say about the "standard" inflationary formula is that it does not 
follow from correct equations. 

It is important to recall that the "standard" inflationary results seem to be an indis- 
pensable tool in cosmology of our days (and possibly in the next Millennium too, see for 
example [41,10]). The expected relative contributions of density perturbations and gravi- 
tational waves to the observed microwave background anisotropics are a subject of active 
study. In the center of discussion are usually the "consistency relations" which state that 
the ratio of density to gravity-wave contributions goes to inflnity when the spectrum of 
perturbations approaches the most favorite, Harrison-Zeldovich form. In reality, as we have 
shown above, these "consistency relations" are simply a maniflstation of inconsistency of 
the "standard" inflationary theory from which they are derived. 

In conclusion, if the "standard" inflationary results are incorrect and cannot be trusted, 
what is the amount of density perturbations that can be generated in the early Universe? 
My part of answer is formulated in Ref. [3,7]. 
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